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Abstract 

Circular Brownian motion models of random matrices were introduced by Dyson and de- 
scribe the parametric eigenparameter correlations of unitary random matrices. For sym- 
metric unitary, self-dual quaternion unitary and an analogue of antisymmetric hermitian 
matrix initial conditions, Brownian dynamics toward the unitary symmetry is analyzed. 
The dynamical correlation functions of arbitrary number of Brownian particles at arbi- 
trary number of times are shown to be written in the forms of quaternion determinants, 
similarly as in the case of hermitian random matrix models. 
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1 Introduction 



Circular ensembles of random matrices were introduced by Dyson as the simplest possible 
models of complex energy spectra[|I], Dyson further introduced a Brownian motion 
model to describe the parametric correlation of the energy levels[§]. His Brownian dy- 
namics is specified by the Fokker-Planck equation 



dp 



(1.1) 

where 61,62, ■■■ ,0^ are the locations of Brownian particles, r is the time variable and 

N 

= -J^logje''^ -e^'|. (1.2) 

j<i 

The stationary distribution of this Brownian dynamics is 

P5 = e-'^^ = n|e'^-e'f (1.3) 

j<i 

satisfying dps/dr = 0. Using a perturbation theory, Dyson showed that it is a natural 
model to describe the parametric eigenparameter correlation of symmetric unitary, unitary 
and self-dual quaternion unitary random matrices, corresponding to f3 = 1,2 and 4, 
respectively. 

In order to solve the Fokker-Planck equation, a transformation of the Fokker-Planck 
operator C into a hermitian operator 7i is useful. The hermitian operator 7i is defined as 

g/3iy/2£e-W2^_l(^_^^^ (1.4) 

P 

with a constant and explicitly written as 

^de]+ 4 ^sin^[(^,-^0/2]- ^ 

It is known as the Calogero- Sutherland Hamiltonian on the unit circle. This Hamiltonian 
has a complete set of orthogonal eigenfunctions 

MOi, ■■■,9^)= e-^'^"Pj^'P\z„ ■■■,zn), z,= e''^ (1.6) 

with eigenvalues -Ek^, i, |H- Here k = {ki,---,^]^) represents a partition with non- 
negative integers ki > K2 > ■ ■ ■ > and P^^/^^(zi, ■ ■ ■ , z^) is a particular polynomial 
known as the Jack polynomial. 

For the imaginary time Schrodinger equation 

dib 1 

= (1.7) 

the Green function solution can be written as 

G(-)(v.i, ...,^j,;e„...,e.;r)=j: Mfhllllf^)^^^ (1.8) 

K \Vk\Vk) 
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where 



J —TV J —TV 



(1.9) 



and ipK is the complex conjugate of ip^- According to (|1.4| ), the Green function solution 
G{ipi, ■ ■ ■ ,ip^;9i, - ■ ■ , On] t) of the Fokker-Planck equation ( |1.1|) is given by 



G(^i 



-(3W{ei,-,eN)/2 



;i.io) 



In this paper we focus on the case /3 = 2 , in which the interaction term in (|1.5| ) 
vanishes, and calculate the dynamical correlation functions for typical initial conditions. 
Then G^^^ is a determinant of one particle Green functions and (|1.10|) reads 



N sm ■ 



j>i sm 



'1.111 



where (z 



w = e^^) 



9(0, 



-(- 

— y 



1/2 oo 

E 



n=— oo 
n 



W 
Z 

--InT 



N even, 
odd 



;i.l2) 



and 



In 



{n - (1/2))V2, N even, 
nV2, A^ odd. 



;i.l3) 



We remark that the following discussion does not depend on the particular form ( |1.13| ) of 
7„, provided they are positive. 

Let us suppose that the initial distribution of the eigenvalues is one of the foUowings 
(note that U{ip) can be a function with a complex value): 



r nf=i u{v,) n 



N 



xn 



Po{ipi, ■ ■ ■ ,V27V) oc < 



assuming A^ even) 



A J < / I II I 

f 1 



;i.i4) 



X 



A^ even, 

S{nN+i)/2])UfJ? I l-e^^^ P, AT odd. 



Here we have introduced a sum J2p over the permutation of (pj^s in order to make the 
initial distribution totally symmetric, [x] is the largest integer not exceeding x. The 
first and second of these initial conditions can be derived as eigenvalue distributions of 
symmetric unitary and self-dual quaternion unitary random matrices, respectively. The 
third one, for which the weight function U{(p) is assumed to be symmetrical about the 
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origin, is an analogue of the eigenvalue distribution of antisymmetric hermitian random 
matrices. 

The probability distribution function p of the eigenvalues is calculated from the initial 
condition and the Green functions as 



■ ■ ■ ; ^A?; Ti; ^1, ■ ■ ■ , ^A?; T"2; ■ • • ; ^f^, ■ ■ ■ , 0^]tm) 

= ^ £ d^? ■ • ■ £ d9%po{9l ■■■,9%)l[ Gi9[-\ • • ■ , 9[, ■■■,9'^;n- 

To = 0. (1.15) 
For even, performing the integration gives 

■ ■ ■ ; ^Ar; Ti; 6*1, ■ ■ ■ , ^'at; T2; ■ ■ ■ ; 9^, ■ ■ ■ , 9'^; tm) 

N N 

^ ^_Af(A^-i)/2 -Q ^-i{N-i)ef/2 J^^Q^ef _ 
i=i 3>i 

M-l 



k=l 



(1.16) 



Here Pf means a Pfaffian. The matrices and are defined as 

^-« = ^(^™ ^«;r„-r„) (1.17) 

and 

i^7" = F(^™ ^«;r„,r„), (1.18) 

where 

F{9,9'-t,t') 

d^' r dypt/(^)t/(^'){^(^,^;^)^(^',¥'V)-^?(^',¥';0^(^,^V)}, 



F{9,9'-T,r 

/TT 
d^? 
-TT 

F(^,^';r,r' 



d^U{vf{g{9, ^; r)^^(^', ^; r') - g{9\ ^- r')-^g{9, v^; r)}. 



dyp[/(^)2— ^{^(^^, r)(7(^', y.; r') - (?(^^', -y;; r')^(e, ^- r)} 



(1.19) 



for each of the three initial conditions ( |1.14| ), respectively. For odd, performing the 
integration similarly yields 

■ ■ ■ 5 ^Af! Tl'^ ^l! ■ ■ ■ ; ^Af! ''"2; ■ • ■ ; ^'f^, • • • , ^jv ; TAf) 

Af AT 
^ ._Ar(Af-l)/2 TT i(Af-l)(?A//2 TT/ ifif _ i9^' 



X Pf 



j=i j>i 

[^^l']j,l=l,2,-,N [fj]j=l,-,N ^ ^ 

-[//]^=i,...,7v 



n det[^?|,-^^ 



jl \j,l=l,2,-,N- 

k=l 



(1.20) 
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Here 

fr = fiOT;rm) (1.21) 

with 

fie;T)= Ui^)gi9,^;T)d^, f{9;r) = UiO)gi9,0;T), (1.22) 

J— TT 

for the first and last of the initial conditions (|1.14|) , respectively. 

Now we are in a position to define multilevel dynamical correlation functions 

P(^i5 ■ ■ • '^mi! ^i; ■ ■ ■ ^Ca! ^2; ■■■■,9^\-- ■ ,9^^j;tm) 
1 {N\)^^ 



X p{9\, ■ ■ ■ ,9lf;Ti;9l,- ■ ■ ,9lf;T2;- ■ 9^^ , ■■ ■ , 9^; Tm) 



;i.23) 



4^ 



where the normalization constant Cat is defined as 

/TT riT PTT PTT 

d9\--- d^jv ■ ■ • / ■■■ 
-TT J —TV J —77 J —77 

X p{9l, ■■■,9\,- n; 9l---, 91- r^; ■ ■ ■ ; <^ r,,). (1.24) 

Let us review the history of the study on multilevel correlation functions of random 
matrices and explain the purpose of this paper. The quaternion determinant formulas 
for the static correlations (M = l,ri = 0) with U{9) = 1 was discovered by Dyson|0]. 
His result was generalized to the case with a real weight function U{9) by Nagao and 
Wadati[P]. The dynamical correlations were first explored for hermitian random matrix 
models by Pandey and Mehta [0, O] and then extended to the unitary matrix case with 



U{9) = 1 by Pandey and Shukla[ll|. They were able to derive quaternion determinant 
forms for equal time (M = 1, < ri < 00) correlation functions. Recently Nagao and 
Forrester succeeded in evaluating the dynamical correlation functions with general M for 



hermitian random matrices [1T2|, [13| . In this paper we deal with unitary random matrices 



and generalize both of Nagao and Wadati's and Pandey and Shukla's results to show 
how the multilevel dynamical correlation functions with general M can be written in 
quaternion determinant forms. 



2 Dynamical Correlation Functions 
2.1 Quaternion Determinant 

Let us begin with an introduction of a quaternion determinant, a determinant of a matrix 
with quaternion elements[0, |T^, |T^. We define a quaternion as a linear combination 

of four basic units {1, ei, 62, 63}: 

g = go + q- e = go + 9iei + ^262 + ^sCs. (2.1) 

Here qo,qi,q2 and are real or complex numbers. The first part qi is called the scalar 
part of q. The quaternion multiplication is associative but in general not commutative: 
the multiplication rule of the four basic units are 

1.1 = 1, l.ej = ej-l = ej, j = 1, 2, 3, (2.2) 



5 



= 62 = eg = 616263 = -1. 
A quaternion q has a dual q defined as 

g = go - q • e. 

A matrix Q with quaternion elements Qji also has a dual matrix Q 
units can be represented as 2 x 2 matrices 

-1 ' 

1 ' 





" -i' 




" i " 


62 


-i 


, 63 ^ 


-i 



(2.5) 



Now wc arc in a position to define a quaternion determinant Tdet. For a self-dual Q 
( i.e., Q — Q ), it is written as 

I 

Tdet Q = El-l)"^"' lii^iabqbc ■ ■ ■ qda)o- (2.6) 
p 1 

Here P denotes any permutation of the indices (1, 2, • • • , N) consisting of / exclusive cycles 
of the form (a — > 6 — > c — > • • • — > d — > a). Note that (—1)^"' is the parity of P. The 
subscript has a meaning that we take the scalar part of the product over each cycle. 
If all the elements of Q are scalars, then everything is commutable and a quaternion 
determinant becomes an ordinary determinant. 

2.2 The Case N Even 

Throughout this paper we adopt a notation that z is the complex conjugate of z (note 
that z* is not necessarily the complex conjugate of z). Let us define {z = e^^, w = e*^) 

K{0;r) = r F{ip,e;T,T)7^R:{ip;T)dip, 

J —IT 

^n{0;T) = r F{ip,9;T,T)V^Rn{^;r)dip, (2.7) 

J— TV 

where R*^{9; r) = + < ^_iZ-''+^ + ■■■ + < and Rn{9; t) = z"" + c„ ^-i^''"^ + • • • + 
c„ are arbitrary polynomials {Rq{9;t) — Ro{9;t) — 1). 

We introduce matrices and 5""^" as 

= D{9]^, 91, r^, r,), /j?" = /(^f , r^, Tn), = ^^^T ^ r^, ^n) (2.8) 

with 

D{9, 9'- T, r') 

(JV/2)-l . 

= E ^^[e-^-'=^^'-^)*:(^;r)*,(^';r') -e-^'=+^(^'-^)*fc(^;r)^^(^';r')], 

(2.9) 



(2.3) 

(2.4) 

[g/j] . The quaternion 



1 
1 



61 
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I{9,e';T,T') 

{N/2)-l . 

+ F{e,9';T,r') (2.10) 



and 



(iV/2)-l , 



Moreover we define 



S^^-g]t^, m>n, 
5]^", m<n. 



(2.11) 
(2.12) 



Then we have the following theorem. 
Theorem 1 

The probability distribution function ( p,.l(j| ) with even can be written as a quaternion 
determinant 

■ ■ ■ ; ^at; ^i; ■ ■ ■ ; ^at; T2; ■ ■ ■ ; 0^ , ■ ■ ■ , 0^;Tm) 



(7V/2)-l 

^-A.(iV-i)/2 -Q ^^.(^^^)Tdet[5n 
i=o 

/i, 1/ = 1, 2, ■ ■ ■ , M. 



:2.13) 



Each block B^'^ is an N x N quaternion matrix. Its quaternion elements are represented 
as 



[2.U) 



Starting from (|1.16|) , we can follow a similar argument as in Refs.[|l^, [T^ to prove Theorem 
1. 

By Schmidt's orthogonalization procedure, we can specify the polynomials Rn{0;T) 
and -R* (^, t) so that they satisfy the following skew orthogonality relation: 



{y/zR*^{9; t), y/wRn{<^; r)) = -{y/zRn{6; r), y/wR*^{i^; r)) = r„(r)(5„ 



where 



{V^R*^i9; r), v^i?:(v^; r)) = 0, {V^Rm{9; r), v^i?„(<^; r)) = 



{f{9),g{^))= d9 d^F{9,^-T,r)f{9)g{^). 



(2.15) 
(2.16) 
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Explicit determinant formula for the skew orthogonal polynomials Rn{6',T) and R*{6,t) 
are given by (n > l)M 



Rn{z) = 



jn n—1 
yU—l jn—1 n—1 



n-1 



n l~n n—1 



jn jn -1 
j'n— 1 J^~^ ~^ 



-1 



and 



r:(z) = v-' 



n J— n—1 —n 



-n+1 



J 



-1 -n 



n Tn—l —n 



J— n—1 J— n—1 1 



J 



jn -n 

n—1 —71 



-n 



J 



-1 



J 



1 1 



J' 



n-1 J"-l 1 



J' 
J' 

J— n—1 n—1 
J^iT- n—1 



J 



-1 n-1 



J' 



n—1 n— 1 



where 
and 



jmn ^ ^_^m+(l/2)^^n+(l/2)^^ 
^ji- 1 n—1 _ _ _ ^n— 1 jn—1 —1 _ _ _ ^n— 1 — n 



n-1 



-n n—1 



JO JO -1 



n 7— n —1 



J' 



-n 



j-n -n 



2.3 The Case N Odd 

Let us next consider the case with odd. We define 

KiO;r) = r Fiip,e;T,T)R:iip;T)dif, 

J — vr 

$n(^;r) = r F{^,e-T,T)Rn{v;T)d^. 

J -IT 

Here Rn{0; r) and -R* (6'; r) are arbitrary functions provided that 

N N 



He 



j(Af-i)ej/2 



X 



-^(Af-l)/2( 



?i;r) R 



(N-l)/2\ 



Rl{e,-T) 
Ro{0i;r) 

Ri{9ur) 



Roi02;r) 
Ri{92;t) 



R(N-l)/2{di;T) R{N -1)/2{G 2] t) 



R] 



{N-l)/2 



{0N;r) 



RUOn^t) 
Ro{9n;t) 
Ri{9N;r) 

R{n-i)/2{Gn]t) 



(2.17) 



(2.18) 



(2.19) 



(2.20) 



(2.21) 



(2.22) 
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Then matrices D'""', /'""■ and S"^^ are introduced as 

= D{9]^, B^- Tm, r.), = I{ef, 91- Tm, r„), = S{ej, e^- r^, r^, (2.23) 

where 



(JV-l)/2 



fc=l ■'^fclTJ 



(2.24) 



/(e,e';r,r') 

(JV-l)/2 



+ 



- E ^^[e~^-^(^'-^^$K^^;r)<ffc(e';r')-e-^'=(^'-^)$.(^^;r)$Ue';r')] 

fc=l ''^fclT J 

<l>o(^; r)f\e'- r') - -^,^9'- r')f\e- r) + F{0, &■ r, r') 



So( 



(2.25) 



and 



{N-l)/2 



+ i?o(e',r')/(e;r). 



As before we set 



nmn 



SJl"", m<n 



(2.26) 



(2.27) 



and find the following theorem. 

Theorem 2 

The probability distribution function (|1.2CI|) with odd can be rewritten as 



■ ■ ■ ; ^at; T-i; ■ ■ ■ ; ^A?; ''"2; ■ ■ ■ ; di^, ■ ■ ■ , ; Tm) 

iN~l)/2 

= r^(^-i)/2so(rM) n r,(rM)Tdet[S'^1, 

/i,i/ = 1,2,---,M. (2.28) 
Each block 5'^'^ is an x quaternion matrix the elements of which are represented as 



(2.29) 



Theorem 2 can be proven from ( l.2(J| ) by proceeding as in Refs.[|T2|, T3 
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We now introduce a set of polynomials J7„(^;r), J7*(^;r) for n > 1 and flo{ 
^11(9; t) satisfying 

{ni{e; r), n„iv- r)) = -{n„i9; r), a;;,(v9; r)) = r^(r)5„„, 

{n*^{9- r), n:(<^; r)) = 0, {QUO; r), Q^^; r)) = 

for m,n > 1 and 

{n:i9; t), noi^; r)) = -{^0(9; r), n^i^; r)) = 0, 

{ni{9- t), fi„(</.; r)) = r), r)) = 

for < n < (A^ — l)/2. Here the bracket is defined in ( p.l6| ). 



(2.30) 



Let us define 



and 



^ f (z'™,fi„((^;r)), n > 0, 
\ (z-,fi:(^;r)), n<0 

(rim(6';r),fi„(v9;r)), m,ra>0, 

(^]„(0;r),fi;(v9;r)), m > 0, n < 0, 

(^];„(e;r),^]„(v9;r)), m<0,n>0, 

(r]:;,(e;r),fi:(^;r)), m, n < 0. 



Then, starting from (2; = e 



i9\ 



ni{9;T) = ai + z, ni{9;T) = al + 
we can recursively construct the polynomials as {n > 2) 



aM^^''\9;r) 



n — 1 
n— 1 



1 

1 

L-1 1 



^n-l -1 



and 



n 



n-lK"i 



T 



J^-n -n+1 



-1 1 



'n— 1 —1 1 



(2.31) 



(2.32) 



(2.33) 



(2.34) 



j^n -n+1 
—n+1 

^1 -n+1 
r —n+1 —n+1 



(2.35) 



J^^iT- n— 1 
rn— 1 n—1 



(2.36) 
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Here 



(n-l). 



-1 

n-l 



^n— 1 n— 1 

^1 n-l 
j^-l n-l 



1 



K' 



-1 



2^-n+l 1 



■n+1 -1 



^0 -n+1 
—n+1 

-"+1 

2^— n+1 —n+1 



(2.37) 



and 



V, 



n-l 



j^n—l n-l 

r —n+1 n-l 



n-l 1 T n—1 —1 



-n+1 1 _^-n+l -1 



—n+1 

^1 -n+1 
^-1 -n+1 



—n+1 —n+1 



(2.38) 



Setting 

no{9;T) = ni^''-'^^'\9;T), (2.39) 

we find that all of the skew orthogonality conditions (|2.3CI|) and (|2.31| ) are satisfied. Note 
that there is an ambiguity in the determination of the skew orthogonal polynomials due 
to the constants a„ and a*. 

Let us now introduce {n > 0) 



SAT 



d9f{9;T)n:{9;T), 
d9f{9-T)nr,{9-T). 



(2.40) 



Then Rn{9; r) and R*{9] r) satisfying (|2]2|) can be constructed as (n = 1, 2, ■ ■ ■ , (A^-l)/2) 



and 



K{9-t) = n*^(9;T)-'^no{9;T), 
R^{9-t) = fi„(0;r)-^(]o(^;r) 



Rl{9-T)=R,{9-T)=n,{9-T). 



(2.41) 
(2.42) 



2.4 Quaternion Determinant Expressions 

Using the orthogonality relations introduced above, it can be readily proven in both the 
cases even and odd that 
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0, 

iDmn 



1 




m < p, p > n, 
m < p, p — n, 





1 



1 



1 



1 



1 



1 



r>mn 



TDmn 



r>mn i r>mn 



TDmn _i_ TDmn 



TDmn _i_ TDmn 



TDmn I TDmn 



m 



p, p>n, 



1 



1 



1 



1 




m — p, p 



1 



1 




TDmn TDmn 



TDmn 




1 



B 



mn 




1 



TDmn 



1 




1 




TDmn 



TDmn 




1 



Moreover we can easily find 



N. 



m > p, p < n, 

m > p, p — n, 

m — p, p < n, 

m > p, p > n, 

m < p, p < n. 

(2.43) 

(2.44) 



We then arrive at the following theorem which has the central importance in the evaluation 
of the dynamical correlations. 

Theorem 3 

In terms of the quaternion elements B^'^, we define quaternion rectangular matrices 



Q 



JL 



TDmn \ 



(2.45) 



V i^T • • • Bf^ , 

Let us further define a self-dual quaternion matrix Q consisting of Qjl 



Q 



/ni nK 



KK 



\ ^nx ni "^nK nx 



Q 



(2.46) 



We denote the j-th row (column) of the m-th row (column) block of the quaternion matrix 
Q as {m,j) row (column). The {m,j) row and {m,j) column contain the variable 9J^. If 
the {rrijj) row and {m^j) column are removed, the resulting smaller matrix Q'^ is still 
self-dual and does not contain the variable 9'^. 
Using the above definitions, we have 

d^^: Tdetg ^{N-n, + l)Tdetg(^ . (2.47) 



/_ 
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The proof of Theorem 3 is found by following the strategy given in Ref . . 

Successive application of Theorem 3 leads to the quaternion determinant expression 
for the multilevel dynamical correlation function as 

p{el ■ ■ • , 9l^-T,- el ...^Bl^-r^-...- , . . . , ^m) = Tdet[i?'^'^K, m,)], 
/i,i/= 1,2,---,M, (2.48) 

where each block B^'^{mfj_, my) is obtained by removing the + 1, + 2, ■ ■ ■ , A^-th rows 
and rriy + 1, rrii, + 2, ■ ■ ■ , iV-th columns from B^^'' . 

3 Dynamical Correlation within the Unitary Sym- 
metry 

In this section we consider the limit Tj — > oo, j = 1,2, ■ ■ ■ , M with all tj — ti, j,l = 
1,2, ■ ■ ■ , M fixed. In this limit it is expected that the dynamical correlations describe the 
Brownian dynamics within the unitary symmetry. We can conveniently take this limit by 
using the following summation formulas. 

3.1 The Case N Even 

Let us express the skew orthogonal polynomials with r = as = e*^) 

n 

j=-n 

R:{9-0) = al^z-^, <„ = !. (3.1) 

j=-n 

Then it can be seen that 

n 

j=-n 
n 

K{9;t) = e-^-"^ E <.z-^e^-^ (3.2) 

j=-n 

with 

r„(r) =e-^-""e-^"+iX(0), (3.3) 
the proof of which comes from an identity 



{V^RUe; r), V^Rniv; r)) = e-^-^e-^-^^^ (V^<(^; 0), V^R^{^; 0)) . (3.4) 
Let us introduce an inverse expansion of p.2|) {n > 1) 



T=0 



n 1 n— 1 

^ng7„+ir = Y,x^,e^^+^^Rj{e-T) + -Yf^nje''-^^R*{0;T), 

j=0 ^ j=0 

n n—1 

,-n^,-.r ^ YK,e'-'^R;iO;T) + zYf^l.e^'-^'^RMr). (3.5) 

j=0 j=0 
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Using ( p.9|) and (|3.5| ), we can readily derive {Xjj = 1, n > 0) 



-7-nT 



^n{e;r) 



27r 



-r„(0) 



u>n 



u>n+l 



2% 



-rn(0) 



u>n+l 



-i^-{l/2)„-7,+ir 



(3.6) 



Putting ( p.5| ) into ( p..l9| ) and comparing the result with ( p.6| ) lead to 



(3.7) 



so that 



I{e,e';r,T') 



OO -| 

k=N/2 ^'^y'^ > 



Substitution of ( |3.6| ) into ( ^.111 ) results in 

S{e, e'- r, r') = 5i(0, 0'; r, r') + 52(^, r, r'), 



where (2; = e , z = e' 



1 /r\V2 



(3.9) 



(3.10) 



n=-(Af/2)+l 



and 



82(6, 6'; t,t') 



27T 



{N/2)-l 00 

E E 



X 



k=0 v=N/2 



(3.11) 



Let us assume that 7^ = 7-n+i and 7„+i > 7^ forn > 1. This assumption is consistent 
with (|1.13| ). In the asymptotic limit Tj 00, j = 1,2, ■• - M with all tj — ti, j,l = 



1, 2, ■ ■ ■ , M fixed, it can be seen from (|3.9|) that 



S{x, y; T, t') ~ Si{x, y; r, r) + ©(e^^^w^'+i^^^/^^"). 



(3.12) 
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We can further derive the asymptotic relations Rn{x;T) ~ 0(1), R^{x;t) ~ 0(1) 
$„(x; r) ~ 0(e-2T"+i^) and $;(x; r) ~ 0(e-2T"+i^) from (|Xg) and (^) so that 



/(x, r, r') ~ 0(e-^('^/2)+i(^+^')), L)(x, y; r, r') ~ 0(e^^/2(^+^')). 



(3.13) 



Therefore in the asymptotic hmit we obtain a quaternion determinant in which S*^" are 
replaced by 

f 5i(^f,0r;^m,r„)-^7r' rn>n, 

and the other elements are set to zero. Then the quaternion determinant becomes an 
ordinary determinant 



cr^^{mi,mi) a^^ {7711,7112) ■■■ cr^*^(mi, mM) 
(T^^(m2,mi) 0-^^(7712,7712) ■■■ a^*'^(?7i2,mM) 

a^^\7nM,mi) ^^^2(^771^,^2) ■■■ a*^*^(mM, ^m) 



(3.15) 



where each block a^^{m^, m^) is obtained by removing the + 1, ?77^ + 2, ■ ■ ■ , A^-th rows 
and rriiy + l, miy + 2, • • • , A'^-th columns from a'^'^. This is the unitary symmetry version of a 



similar determinant formula [|1^, ^ known for the Brownian dynamics within hermitian 
symmetry. 



3.2 The Case N Odd 

In the case with odd N, we introduce expansions (77 > 1, 2; = e*^) 



j=-n+l 
n 



j = -71+l 



with = ttnn = 1 and 



(N-l)/2 

fio(^;r) = e-^°^ ^ ao.zV^". 

j=-{N-l)/2 



As before it can be easily proven that 



r„ r = e 



-^-""e-^""r„(0) 



and 



s„(r) = e-^-.„(0), <(r) = e-^--<(0). 



(3.16) 



(3.17) 



(3.18) 
(3.19) 
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Let us now define 

n:(^;r) 



— < 



0, 



< n < (AT - l)/2, 



27r 



(0;((^;0),Oo(<^';0))|^^oTr(e;T), n> (iV + l)/2, 



— < 



< n < (iV- l)/2, 



-7nT 



I 27r 



where 



(n„((^;0),no(<^';0))|^^oT(^;T), n> (7V + l)/2, 



1=1 



The inverse expansions (n > 1) 



"A 



.ng7nr ^ ^A„,e^^^^(^;T) + ^//„,-e^--0*(^;T), 



i=o 

"A 



j=0 j=l 



and 



where 



(Ar-l)/2 (-'V-i)/2 

^ Ao,e^^-"Jl,(^;T)+ 5^ //o,e^-"ll*(^;^)' 

j=0 j=l 



nx = + 1 = n, if n > (A^ - l)/2, 

nx = n^ = (7V-l)/2, if n < (TV - l)/2, 



lead to 



= < 



-7-„r!n(0) 

27r 

-7-.r!n(0) 

27r 



n„(^;r) 

7.r^n(0) 



(3.20) 
(3.21) 



(3.22) 



(3.23) 



(3.24) 



Ll=0 1=1 



l=n 



l=n+l 



l<n< {N-l)/2, 
n> {N+l)/2, 



— e 



— e 



27r 

-7„r^n(0) 

27r 



E 



^ ^-'e-^'>„ + ^^'e-^-'-AL 

=n+l l=n 



l<n< (TV - l)/2, 
n > (iV + l)/2. 



(3.25) 
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Using the formula ( p.25|) for n*(^; r) and n„(^; r), we can readily derive 

OO . 1 

/(^' ^) = E rTTvKl^; r).„(r) - n„(0; r)<(r)] + -.o(0)T(^^; r). 

n=l ' / 

Substituting and (jX^) into ( CT) yields 

8(9, 9'; T, t') = 3,(9, 9'- r, r') + S-,{9, 9'- r, r'), 
where [z = e*^, = e*^') 



Sii9,9';T,T') 



27r 



{Ar-l)/2 

E 

n=-(Af-l)/2 



(3.26) 



(3.27) 



(3.28) 



and 



1 



(7V-l)/2 



X 



S2i9,9';ry) = - E E 

^" fc=0 i/=(Ar+l)/2 

(z-'^e-^'-^A^fc + ;2'^e-^-Xfc)e^'="'^A.(^'; r') 

e^--'fi*(0';rO 



+ 



vk 



+ 



So(^ 



n^{9':T') ~ 



(iV-l)/2 



J\^{9\t)'-"^^^ 1 



fc=i 



'^fc(^' 



^o(r') 



E 



1 



fc=(Af+l)/2 



K(^;r)^/c(r)-n,(0;r)4(r)] 



(3.29) 



where we set q = fi* q = 0. Finally we put and ( ^^231) into (ITIgl) to find 

Ii9,9';T,T') 



E 



+ 



fc=(7V+l)/2 

no(^;r) - 



rfc(rO 

oo 

E 



e-7-.{r'-r)n*(^. ^)n,(^^'; r') - e-^'=(-'--)nfe(^; T)Uli9'; r'] 
1 



fc=(Af+l)/2 
oo 

E 



So{r 



^o(r') fc=(^i)/2 ^fc(^) 



-[n*(0';r')^;.(T)-n.(0';r')4(T')] 



nn(^';r') 1 



[Uli9;T)sk{T)~U,i9;T)sliT)]. 



(3.30) 



Assuming that 7„ = 7_n and 7n+i > 7n for ^ > 0, we can now take the asymptotic limit 
Tj —>■ oo, j = 1,2, ■ ■ ■ M with all tj — ti, j,l = 1,2, ■ ■ ■ , M fixed. Let us first note that 
n„(x;r) ~ 0(1), fi;(x;r) ~ 0(1), n„(a;;r) ~ 0(e-2^"") and n;(x;r) ~ 0(e-2^"") for 
n > 1. We can further obtain an estimation 

Qoi9; r) ~ 0(e-(^o-^(^-i)/^)"), no(^; r) ~ 0(e-(^°+^(^+i)/2)-). (3.31) 

Then it is straightforward to find 

S{x, y; T, t') ~ Si{x, y; r, r') + 0(e-(^('^+i)/2-T('^-i)/2)"), (3.32) 
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I{x, y; r, r') ~ 0(e-'^(^+i)/2 D{x, y; r, r') ~ 0(e'^(^-i)/2 (3.33) 

Therefore and can be set to zero in the asymptotic hmit and, as in the case N 
even, a determinant expression ( |3.15| ) results. Here the matrix elements o"^" are again 
defined by ( p.l4| ), although the definitions of 5*1 (6', 9'; r, r') and g{9, ip; r) are different from 
those in the case N even. 



4 The Case U{0) = 1 

In this section we deal with the simplest case U{6) = 1 and explicitly give the formulas 
for skew orthogonal polynomials. Though this case was already treated by Pandey and 
Shukla, our result is more general than theirs because they evaluated only the equal time 
correlations. 



4.1 The Case N Even 

(1) Symmetric unitary initial condition 

For the bracket defined in (|2.16|) corresponding to the first of the three initial conditions 
( |1.14|) , we can easily derive {z = e*^) 



^^m+(l/2)^^n+(l/2)^ 



=0 n+(l/2) 



which yields 



Rnie;0) = z^, R:i9;0) = z-^, r„(0) 



n+{l/2)' 



(2) Self-dual quaternion unitary initial condition 



Let us next consider the second of the initial conditions (|1.14|) . Noting 



[z 



m+(l/2) ^n+{l/2)\ 



r=0 



we obtain 



Rn{e;0) = z^, R:{e;0) = z-^, rn{0) = Am (n + 

(3) An analogue of antisymmetric hermitian initial condition 
The third of the initial conditions ( |1.14| ) yields 



^^m+(l/2)^^n+{l/2)^ 



T=0 



' Tii, if n — m is odd and positive, 
0, if — m is even, 
— vri, if n — m is odd and negative. 



which implies {n > 1) 

i?„(0;O) = ;2" + z-", i?:(^;0) = 2-" + z", r„(0) = 27r^ 

and 

i?o(^;0) = l, ro(0)=7ri. 



(4.1) 



(4.2) 



(4.3) 
(4.4) 



(4.5) 

(4.6) 
(4.7) 
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4.2 The Case N Odd 

(1) Symmetric unitary initial condition 

In N odd case of the symmetric unitary initial condition, we obtain 

4-7r? 

(^'",'2")Uo = — ^-+no, m,n^O (4.8) 
n 

and 

4-7r? 

(Ai)u = — (-ir- (4.9) 

Then, putting * = 0, we can derive (n > 1) 

A.'jrj 

O„(e;0) = ^^ (]:(^;0) = z-", r„(0) = — (4.10) 

and 

(JV-l)/2 

^o(^;0)= E (-^)^'. (4.11) 

j=-(Ar-l)/2 

(2) An analogue of antisymmetric hermitian initial condition 

An analogue of the antisymmetric hermitian initial condition is handled by considering 
the inner products 

{ni, if n — m is odd and positive, 
0, if n-m is even, (4.12) 
—Tri, if n — m is odd and negative. 

We fix the constants as a„ = a* = (—1)"+^ and find {n > 1) 

Qn{0;0) = z"" + z-''^^ Ql{e;0) = z-'' + z''-^ r„(0)=27rz (4.13) 

and 

n,(9; 0) = (_i)(^-i)/2 1 (^(iv-i)/2 ^ ^-(iv-i)/2^_ 

5 Conclusion 

In this paper we have derived quaternion determinant expressions for the dynamical mul- 
tilevel correlation functions for Dyson's Brownian motion of eigenparameters (energy 
levels) toward unitary symmetry. As the initial conditions, we assume eigenparameter 
distributions of symmetric unitary, self-dual quaternion unitary and an analogue of the 
antisymmetric hermitian random matrices. Any of the quaternion elements is represented 
in terms of four functions D{9, 6'; r, r'), I{6, 6'; r, r'), S{6, 6'; r, r') and g{6, 6'; r — r') ap- 
pearing in the two level dynamical correlation functions. Therefore analyzing the two 
level dynamical correlations is enough to understand the behavior of all the multilevel 
correlations. The two level dynamical correlations are investigated in Refs. |p]|, pi] , ^ 
and some asymptotic results in the limit N oo are already known. 
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